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The dynamics of weakly coupled, non-abelian gauge fields at high temperature is 
non-perturbative if the characteristic momentum scale is of order |k| ~ g'^T. Such a 
situation is typical for the processes of electroweak baryon number violation in the 
early Universe. Bodeker has derived an effective theory that describes the dynamics 
of the soft field modes by means of a Langevin equation. This effective theory 
has been used for lattice calculations so far [18, 19]. In this work we provide a 
complementary, more analytic approach based on Dyson-Schwinger equations. 
Using methods known from stochastic quantisation, we recast Bodeker's LANGEVIN 
equation in the form of a field theoretic path integral. We introduce gauge ghosts in 
order to help control possible gauge artefacts that might appear after truncation, and 
which leads to a BRST symmetric formulation and to corresponding Ward identities. 
A second set of Ward identities, refiecting the origin of the theory in a stochastic 
differential equation, is also obtained. Finally Dyson-Schwinger equations are 
derived. 



I. INTRODUCTION 

Lattice calculations are an ideal tool to extract with a minimum of theoretical prejudice 
a specific piece of information from a given theory. However, in a sense, they are kind of a 
'black box' that give the answer but hide the way how that answer comes about. 
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The aim of this work is to provide a complementary, more analytic approach to the non- 
perturbative physics encoded in Bodeker's effective theory [1]. The emphasis thereby lays 
not primarily on the accuracy of the results where it is hardly possible to beat the lattice 
calculations. Our aim is to provide a tool for a deeper understanding of what is really going 
on in the non-pcrturbativc sector of hot non-abelian gauge theory and during creation of 
baryon number. In particular, it allows for an analytic study of the sphaleron rate [2, 3] 

y— +00 1— »oo vt 

Such a deeper understanding is not only important for baryogcncsis. Magnetic screening and 
the corresponding identification of a magnetic mass are of quite general theoretical interest 
with applications also in the field of the quark-gluon plasma [4-7] . 

We base our analysis on Bodeker's effective theory, despite the fact that Bodeker has 
also derived a generalised Boltzmann-L ANGEVIN equation which is valid to all orders in 
[log(l/5f)]~^ [9], and of which Bodeker's effective theory is merely the leading logarithmic 
approximation. We choose this approximation because the more general Boltzmann— 
Langevin equation not only is far more complicated, but is also not renormalisable by 
power counting [10]. The effective theory on the other hand is ultraviolet finite, and is 
known to still be valid at next-to-leading logarithmic order provided one uses the next-to- 
leading logarithmic order colour conductivity a [11]. 

The key idea of this work is rather simple. Bodeker's effective theory has the form of 
a Langevin equation. It is well-known from stochastic quantisation that a Langevin 
equation can be recast in the form of a path integral [15-17]. This path integral then 
can be reinterpreted as the functional integral formulation of an Euclidean quantum field 
theory with some 'strange' action. In this way, one gains access to all the powerful methods 
developed in QFT. Specifically, it is possible to derive the Dyson-Schwinger equations 
of the theory, offering an approach to the non-perturbative sector that is independent from, 
and complementary to, the existing lattice studies. 

On the way to this goal a couple of obstacles have to be overcome. These are mostly 
related to the peculiar role played by gauge invariance in the context of stochastic quanti- 
sation and Bodeker's effective theory. A thorough understanding of this role proves to be 
essential in pursuing our aim. 

The outhne of this work is as follows. Section II is devoted to the transcription of 
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Bodeker's theory in path integral form. Prom this path integral one could proceed to de- 
rive the Dyson-Schwinger equations, and in principle, could gain access to the non- 
perturbative sector of the theory. At the end of the day however, one will be forced to rely 
on a certain truncation scheme to extract any concrete results from the equations. This 
truncation may introduce a possible gauge dependence and thus may render the results 
worthless. To keep control over the gauge dependence, it is therefore necessary to generalise 
Bodeker's equation from Aq — gauge to a more general class of gauges before applying the 
formalism. 

Gauge fixing in a stochastic differential equation is quite delicate. One has to make sure 
not to destroy the Markovian nature of the equation. Applying methods developed in 
the context of stochastic quantisation [12], we introduce a gauge fixing term into Bodeker's 
equation thereby achieving the desired upgrade to a general class of flow gauges. 

In Section III we argue that any physically reasonable truncation of the Dyson- 
SCHWINGER equations requires the introduction of gauge ghosts. In the full, untruncated 
theory gauge ghosts are not necessary, which is generally true in stochastic quantisation 
[12, 16]. As was shown in Ref. [12], gauge ghosts can be introduced in stochastic quantisa- 
tion in order to establish a gauge BRST symmetric formulation. 

We carry out this program in the case of Bodeker's theory and derive the Ward- 
Takahashi identities corresponding to the gauge BRST symmetry of the action. These 
should be respected by the truncations to be used. 

The gauge Ward identities are not the only restrictions to be observed. A second class of 
Ward identities exist, that are related to the characteristic structure of the theory reflecting 
its origin in a stochastic differential equation. This characteristic structure can as well be 
expressed in the form of a BRST symmetry by introducing another kind of ghost flelds, 
referred to as equation of motion (EOM) ghosts in this work. Introducing the gauge ghosts, 
however, destroys this second stochastic BRST symmetry. Nevertheless, it does not change 
the physical contents of the theory. The stochastic BRST symmetry is only an elegant 
way to express this structure. By directly referring to the underlying physics, it is still 
possible to derive the corresponding stochastic Ward identities. They provide a second set 
of restrictions to be imposed on the truncations. 

In Section IV we derive the Dyson-Schwinger equations of Bodeker's effective theory. 
In combination with the gauge and stochastic Ward identities of Section III, this constitutes 
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an independent approach to the non-perturbative dynamics of the soft, non-abehan gauge 
fields encoded in Bodeker's effective theory. 

In Section V we summarise and discuss our results. Appendix A shows the explicit calcu- 
lation of some of the Jacobians encountered in this work. We have included this Appendix 
in order to make our presentation more self-contained. The Feynman rules corresponding 
to our field theoretic transcription of Bodeker's effective theory are listed in Appendix B. 
Finally, in Appendix C, we present explicit identities for the lower n-point function following 
from the general Ward identities. 

II. PATH INTEGRAL FORMULATION OF BODEKER'S THEORY 

A. Transcription to a Path Integral in = Gauge 

According to Bodeker's effective theory the dynamics of the soft modes of the gauge field 
is described to leading logarithmic order by the Langevin equation [1] 

D"'' X B'' -h ak^ = C (2) 

which is written in = gauge and where ^ is a gaussian white noise stochastic force. 
The stochastic force field incorporates the infiuence of higher momentum modes and has the 
correlator 

(C(i,x)C''^(i',x')> = 2aTS'^5''^5{t - t') 5^-^(x - x') (3) 

reflecting its gaussian white noise character. Here and in the following, the number of 
spacial dimensions is D — 1 = 3, however, we leave D unspecified to allow for dimensional 
regular isat ion later. The only physical parameters entering Eqs. (2) and (3), and therefore 
the effective theory, are the temperature T, the colour conductivity cr, and the self coupling 
of the gauge field hidden in the definition of the covariant derivative D"^ = 5"'^V — gf'^'^A^. 

The procedure of reformulating a Langevin equation like Eq. (2) in the form of a field 
theoretic path integral is well-known [15-17]: According to Eq. (2), the gauge field evolves, 
starting from certain initial conditions, under the influence of the stochastic force. An 
arbitrary observable of the theory then is deflned by some functional of the gauge fleld F[Pl\ 
and given by the expectation value of that functional with respect to the possible realisations 
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of the stochastic force 

{F[A])=Jl)CF[A%MC]=J'DCF[A%]]expJ^-^Jdtd^-'xC{t,^ (4) 

Here we have denoted by A^[C] the solution of Eq. (2) for a specific choice of the stochastic 
force and the given initial conditions. 

To proceed and recast the effective theory of the gauge field into a form resembling the 
path integral formulation of an 'ordinary' quantum field theory we would rather like to have 
a path integral running over the gauge field than running over the stochastic force. This can 
be achieved by inserting unity in an appropriate way. In fact, one has 



J d{E - C) = J DA Det (^f^) 5(E[A] - C) 



(5) 



(i.c.) 

where we choose the functional E[A] as the left-hand side of Eq. (2) 

E"[A] = D"^x B'' + (tA" (6) 

The invertibihty of E[A] is essential to justify the change of variables in Eq. (5). It follows 
from the parabolic nature of the expression and from the restriction to those gauge field 
configurations in the second path integral satisfying the initial conditions. 

Because Eq. (5) holds independently of it can be inserted into the path integral (4). 
The delta function then assures that only those gauge field configurations contribute to the 
integral that obey E[A] = Due to our choice of E[A], however, this is identical to the 
condition A = A^[C]. Thus, after inserting the delta function we may replace A^[^] in the 
path integral simply by the integration variable A, and we are left with 

(F[A]) =y 2?C ^[C] / m Det (^^^ 5(E[A] - C) F[A] (7) 

(i.c.) 

Moreover, the restriction to field configurations obeying a specific set of initial conditions can 
be dropped if these initial conditions are specified at i = — oo. This is a consequence of their 
transversal component always being damped and the fact that any longitudinal contribution 
drops out whenever a gauge invariant observable is calculated. In case of a gauge variant 
quantity, however, a damping of the longitudinal component can be achieved by introducing 
an additional gauge fixing term into the Langevin equation [16]. This will be necessary 
anyway in the following section in order to generalise from Aq — gauge. Henceforth, we 
will therefore drop the restriction on the path integration in Eq. (7). 



6 



At this point, one has two choices. One possibihty is to proceed by doing the C integral 
with the help of the delta function. This results in a theory containing only the gauge field 
(and perhaps some additional ghost fields to be introduced later), however, at the expense of 
rather complicated interactions: the functional E[A] shows up as argument of the gaussian 
probability distribution, and since E[A] contains terms up to A^, the action would inherit 
vertices of up to sixth order. 

To avoid this situation, we instead choose to introduce an additional auxiliary field A to 
represent the delta function 

.(E[A,-C)^/mexp A»,E«[A,-0} (8) 

In this way, one can still perform the C integral that becomes gaussian, thereby eliminating 
the stochastic force field from the theory. One obtains 

(F[A]) = J DADA Det (^^|^) ^[A] e"^!^'^] (9) 

with 

S[A, A] = Jdx [(jT A"-A" - iA'^- E"[A]] (10) 

The determinant in Eq. (9) need not be taken into account since it can be shown to be 
a constant in dimensional regularisation (see Appendix A for an explicit calculation). We 
could, nevertheless, introduce a ghost representation of the determinant referring to the 
corresponding ghost fields as equation of motion (EOM) ghosts in the following. As a 
benefit of doing so the action (10) would be endowed with a BRST symmetry, allowing 
to easily obtain a kind of Ward identities (so-called stochastic Ward identities) reflecting 
the origin of the theory in a stochastic differential equation. Since it is desirable to obtain 
as many non-perturbative identities as possible in order to find a judicious ansatz for the 
truncation of the Dyson-Schwinger equations, introducing EOM ghosts, at first, seems 
the natural way to proceed. 

However there is another type of Ward identities related to gauge invariance. Unfortu- 
nately, the gauge ghosts to be introduced to obtain these gauge Ward identities will break 
the stochastic BRST symmetries. So, instead of introducing EOM ghosts now, we will later 
introduce gauge ghosts in order to obtain the gauge Ward identities. The stochastic Ward 
identities will be derived without the help of a BRST symmetry by directly referring to 
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the fundamental structure of the theory that reflects its origin in a stochastic differential 
equation. 

For now, absorbing the constant determinant in the measure, we are left with 

(F[A]) mmF[A] e-^[^'^l (11) 
where the action S is given by Eq. (10). 



B. Upgrading to k Gauge 

Bodeker's theory is written in Aq = gauge, and so is our transcription as field theoretic 
path integral so far. At the end of the day, however, we will be forced to use an approximation 
to solve the non-perturbative equations obtained, e.g. Dyson-Schwinger equations, and 
this approximation might introduce gauge artefacts into the calculation. In order to allow 
some control over the gauge dependence of the results, we need to base our derivations on 
a reformulation of Bodeker's equation in a more general gauge. 

In [12], Zinn- Justin and Zwanziger have shown that adding a term to Eq. (2) that is 
tangent to the gauge orbit 

X + (7(A" + T>'^v^[A]) = C (12) 

has no effect on expectation values of gauge- invariant objects of the form F[A]. This is not 
the most general modification of Eq. (2) which leaves expectation values of gauge invariant 
objects unchanged [14], but it suffices for our purposes. As long as f"[A] contains no time 
derivatives, the added term has no effect in calculations of gauge invariant objects. 

We can reformulate this fact in a different way: Since the non-abelian electric field is 
given by E" = —A" — JD^^A^^, one may rewrite Bodeker's equation in the compact form 

D"^ X - (jE'^ = C" (13) 

which then may be interpreted in any of the so-called flow gauges A"'^ = ^"[A] with no time 
derivatives allowed inside the functional v"[A]. 

The restriction that v"'[A] does not contain time derivatives plays a more substantial role 
in our context than in the context of stochastic quantisation which was the object of Zinn- 
Justin and Zwanziger: In stochastic quantisation the time variable describes a fictitious time 
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that is introduced only as a device to reinterpret a given Euchdean quantum field theory as 
the limit of a stochastic process for large values of the fictitious time [16]. Absence of time 
derivatives in stochastic quantisation therefore means absence of derivatives with respect to 
fictitious time and does not pose any restrictions to usual time derivatives. In our context, 
on the contrary, time is the real, physical time and the restrictions above narrow down the 
class of possible gauges leading to a well defined Langevin equation. 

Moreover, because of the different role of the time variable, we also have a component of 
the gauge field that is associated with the t variable of the Langevin equation. In stochastic 
quantisation this is not the case because t is fictitious and the time associated with Aq is 
just the zero component of the Euclidean x vector. To cope with this different structure, to 
some extent will demand a generalisation of the proof of Zinn- Justin and Zwanziger. 

In effect, we not only have to prove that gauge invariant objects of the form F[A] are 
left invariant by the introduction of the term i'"[A], as was shown in [12]. Instead we have 
to prove the following: Given Bodeker's equation in the form (13) and a gauge invariant 
functional ^[^4°, A], then any choice of a fiow gauge leads to the same result. Or put in 
different words, calculating (F[v[A], A]) by means of the equation Eq. (12) gives always 
the same value, independent of ^[A]. 

We now proceed in a similar manner to [12]. Let us consider the left-hand side of Eq. (12) 
where we add a small variation of the v"[A] term. We evaluate this expression for a gauge 
field that is subject to an arbitrary, infinitesimal gauge transformation A'" = A" + D^^'o;^ 
and find 



j3/ab ^ g/fe ^ ^^^/a ^ D'ab^b[A'] + D'"''5v^[A']) (14) 

= {S'''' + gr''^u;'') D^'^x B'^ + (7(A^ + D^V[A])| + aB"' 



+ [H[A]u;f + Sv'[A] 



dt 



Here we have used 

A'" = {6"'' + gr^'uj") + D"^ — (16) 

i.e. the product D"'' x transforms covariantly whereas the transformation of A" has a 
covariant and non-covariant contribution. In the same way we have split the transformation 
of v"'[A\ into a covariant and non-covariant part: Starting from 

v^[A']{t,^) = v^[A]{t,^) + jd^-hj^-^^^^^SA^M (17) 
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we have indeed 

v"[A'] {t, x) = ((5"'' + gr^^Lu") v^[A] {t, x) + [H[A]iu]''{t, x) (18) 

where SA^^ = D^'^lu'^ has been used and we have introduced the abbreviation 

[i/[A]a;]"(t,x) = /^"""V^^^^^ (Dfa;^)(t,y) -5rS^[A](t,x)u;'^(t,x) (19) 

Note that the functional derivatives in Eqs. (17) and (19) are only with respect to a spacial 
variation because f"[A] docs not contain any time derivatives (otherwise wc would also have 
to integrate over time). Let us give the explicit form of this somewhat frightening expression 
for H[A]uj in the case of the choice v^[A] — — • A". One simply obtains 

[H[A]uj]%t, x) = -- (D"'' • Vu'){t, x) (20) 

K 

Finally Eq. (18) leads to 

D'"^^[A'] - {S''^ + gf^^u") Ti^%'^[A] + D"'' [H[A]ujf (21) 
where it was used that uj is infinitesimal and of course 

D'"^5v''[A'] = D"^(5v^[A] (22) 

because 5v is infinitesimal itself. 

Let us now come back to Eq. (14) and its meaning. Suppose the gauge field, before the 
gauge transformation has been performed, was a solution of Bodeker's equation with the 
v°'[A\ term present, but without the additional 5t;"[A] term. In other words, the original 
gauge field was a solution of Eq. (12). We can then replace the first square bracket on the 
right-hand side of Eq. (14) by the stochastic force and find 



+ [H[A]u;f + dv'[A] 



(23) 



dt 

This means, if we subject the original gauge field to an arbitrary, infinitesimal gauge trans- 
formation with parameter cu, then the gauge transformed field will be a solution of Eq. (23), 
i.e. of the original equation with v replaced hy v + Sv and the stochastic force transformed 
in the same way as the gauge field . . . but with an ugly additional term on the right-hand 
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side. However, one can play a dirty trick: What was said so far was true for an arbitrary 
gauge transformation. But if we demand a; to be a solution of 

a, ,b 

^ + [H[A]u;]' + Sv'[A] = (24) 
then the square bracket on the right of Eq. (23) will vanish and we finally arrive at 

Q/ab X g/b ^ ^(^/a ^ D'ab^b^A'] + B"'^6v^[A']) = C (25) 

However, there is a certain subtlety that we want to draw attention to. To clarify this point, 
let us once again repeat the line of reasoning: Starting with a gauge field being solution of 

D"*" X B*" + (7(A" + D«^^[A]) = C" (26) 

we search for a gauge transformation cu that obeys 

"-^ + [H[K]ujf + 5v'^[A]^Q (27) 

(and we can always find such an uj because (27) is a linear, inhomogeneous equation with 
given inhomogeneity 5i'"[A]). Then the gauge field transformed with this uj, A'" = A" + 
D"^a;^, is a solution of the original equation with v replaced hy v -\- 5v and the stochastic 
force also transformed by the same uj 

Q/ab ^ g/b ^ ^^j^a ^ + Ti'^^Sv^A!]) = C'" (28) 

The subtle point is the following: The original gauge field A is a solution of Eq. (26) and thus 
depends on the stochastic force C, of course. But A is an input of Eq. (27) that determines 
UJ. Therefore, uj via A too depends on C,- As a consequence of this, C,' inherits a non-trivial 

dependence on C,: The stochastic force C,' not only depends on C, because it is the gauge 
transform of C, but also because the gauge transformation itself depends on C, 

C^{5^' + gr'''oj%])C (29) 

We denote by A^[^, t;, Aini] the solution of Eq. (12) for the specific reahsation ^ of the 
stochastic force term and initial conditions Aini. Correspondingly, let A*[^,i> + 5v,Ai-a^ 
denote the solution of this equation with v replaced hy v + 5v and for the same stochastic 
force and initial conditions. We can then express the contents of Eq. (28) in this new notation 

Arc V + 5v, -Ai„i] = -A^[C, V, Aini] (30) 
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where the superscript u indicates gauge transformation with the special parameter uj corre- 
sponding to the solution on the right-hand side via Eq. (27). 

After these preparations we can now show that gauge invariant expectation values {F[A^, A]) 
are independent of the choice of f"[A]. To this end, let us write the gauge invariant observ- 
able as functional of the non-abelian electric and magnetic field 

(31) 

B« = V X A" + i gf^-^^A^ X A^ 

We then have 
with 

K+5v[M = -A" - D"^^[A] - D«^5v''[A] (33) 
and B^_|_5^[A] = B^[A] as in Eq. (31). Changing variables according to Eq. (29), one obtains 

(F[E,B]),^,, = j 2)CDet (^^j e[X] F[E„+,,[A], B,+,,[A] (34) 

We now use independence on the initial conditions, the transformation property (30), gauge 
invariance of q[C,] and finally the fact that the determinant is unity (shown in Appendix A). 
This all together leads to 

(F[E,B])^^,„ = ^[C] F[^v+5.rM.'^v+5.m]^=^s^^^,,^.^^^ (35) 

Taking into account the transformation properties (16), (21) and (22), we find 



K^SvW = (-E,[A])'^-D«''[^|^ + [i^[A]^]" + 5^''[A]J =(-E,[A]) 
K^SvW = ("B,[A])« 

and thus 



(36) 



(^[E,B])^^,^ = j^C q[C] ^["E4A],-B,[A]]^^^,j^^^_^^^^j = (F[E,B])^ (37) 
because F[E, B] is a gauge invariant functional. 

Consequently, we have shown that Bodeker's equation in = gauge 



D^^'x B^ + (7A« = C" (38) 
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can equivalently be formulated in any flow gauge 



D"*" X B'' + (7(A" + D"^''[A]) = C 



■a 



(39) 



without any time derivatives allowed inside the functional t''*[A]. We will henceforth use the 



for v"[A], since it has the lowest order in A, preserves colour invariance, and with k > the 
term D"^i>^[A] provides a globally restoring force along gauge orbits [13], while at the same 
time having the correct dimensions. 

III. BRST SYMMETRIC ACTION AND WARD-TAKAHASHI IDENTITIES 

We have argued that in order to derive any reliable statements from our theory, it is 
essential to gain some control over the gauge dependence possibly introduced by the trun- 
cation of the Dyson-Schwinger equations. This was our main motivation to generalise 
Bodeker's equation from = gauge to a more general class of flow gauges. In addition 
to this, the corresponding introduction of a gauge-fixing force has a welcome side-effect: It 
solves at the same time the problem of undamped longitudinal components of the initial 
gauge field configuration. 

However, the detection of an unphysical gauge dependence is not what we really want; in 
fact, we would rather like to avoid it. The ultimate goal is to construct a truncation scheme 
that is physically reasonable and does not (or, reahstically speaking, only slightly) violate 
the gauge symmetry. 

To this end, we need identities expressing the gauge symmetry on the level of n-point 
functions, i.e. we need the Ward-Takahashi identities of the theory.^ 

Any physically reasonable truncation will have to respect these identities. Besides this 
conceptual importance, we may also hope that some of the Ward identities to be derived 
in the following will be of some practical use in solving the Dyson-Schwinger equations: 
In ordinary QCD, for instance, the full gluon propagator in covariant gauge is restricted to 

^ In the non-abelian context, these identities are often referred to as Slavnov- Taylor identities. However, 
following the terminology of Ref. [12], we denote these identities as gauge Ward identities in analogy to 
the stochastic Ward identities also encountered in this work. 



special choice ^4"° = v"[A] 



= — • A" and refer to it as k gauge. This is a natural choice 
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being purely transversal as a consequence of the Ward identities. This leads, of course, to a 
great simplification in the Dyson-Schwinger equations of QCD. 

In this section, we study three different kinds of non-perturbative identities: gauge Ward 
identities, i.e. Slavnov— TAYLOR identities; stochastic Ward identities; and ghost number 
conservation. 



A. Constructing a BRST Symmetric Action 

In Section II B, we saw that Eq. (12) transforms covariantly only under a restricted class 
of gauge transformations. Obtaining the gauge Ward identities with this restriction turns out 
to be rather cumbersome. Instead, we will raise to life the gauge parameter lo by introducing 
into the theory an additional (Grassmann valued) field that realises the constraint on the 
gauge transformations. The resulting action will be endowed with a BRST symmetry, and 
we will be able to obtain the gauge Ward identities in a straight-forward manner. 

Setting (^I'^fA] to zero in Eq. (24) we see that Eq. (12) transforms covariantly under gauge 
transformations which obey 

^ + [H[A]u;f = (40) 

Note that the introduction of v°'[A] does not restrict the gauge group any further than it 
already would be. Even without the extra term, the gauge transformations would have to 
be restricted in order for Eq. (12) to be gauge covariant. 

The restriction in Eq. (40) can be taken into account in the path integral in the following 
way. Define a term 7"[u;, A] from the left-hand side of Eq. (40), which for our choice of the 
functional ^"[A] takes the form 

a, .a 1 

7"[a;,A] = ^--D'^''-Va;'' (41) 

Perform a change of variables from 7 to a; in the following Grassmann integral representation 
of unity 

j ' j Det(^) ^ dt ^ ^ ^ ' 

Since the determinant is Grassmann even, it no longer depends on uo and it can be pulled out 
of the integral. The determinant is a constant, and can be calculated in a similar manner 
to the determinant in Eq. (9); see Appendix A for an explicit calculation. 
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Inserting the integral representation of the Grassmann delta function 

5(7) = y^c^ exp i^Jdx LU^ix) 7"(a;)| (43) 

and absorbing the constant determinant into the measure, we find the identity 

1 = J DuDu exp I Jdx (^5"''^ - ^ D"''- \^^uj\x)^ (44) 

which holds independently of the gauge field A. Therefore, it can be inserted into the 
analogous of the path integral representation of the generating functional, Eq. (11), based 
on the generalised version of Bodeker's equation (12). This leads to 

Z[J] = J TlAVXDuj'Duj exp|-5[A, A,cu,cD] + jdx J''(a;)A"(a;)| (45) 

with the action now given by 

S[A,\uj,uj\ =5(°)[A,A] + 5(°°)[A,u;,a;] (46) 

where 5**°^ [A, A] is the generalised contribution of the dynamical fields 

5'°)[A,A] = jdx cTTA"-A"-iA"-(D"^xB'' + (7(A«-iD"''V-A^)) (47) 

and 

5(°°'[A,a;,a;] = ^c^a; [- a;"cO» + ^ tD"D'^^ • Vo;''] (48) 
is the new contribution containing the gauge ghosts u) and uj. 

B. Gauge Ward Identities 

The Slavnov-Taylor identities can be derived by noting that the action (46) is invariant 
under the following BRST transformation 

5eA^{x) = T>''\x) euj^{x) Se^'^ix) = \gf''^^euj\x)uj\x) 

(49) 

5^X{x) = gf"-^''euj''{x))^{x) S^ui'^ix) = gf''^''£uj%x)uj\x) + ieaTi''^{x)-)^{x) 

where £ is a constant Grassmann parameter. It is convenient to introduce the finite BRST 

operator s such that the result of s acting on a functional of the fields A, A, u and uj is 
defined as (left) derivative with respect to the parameter e of the variations in Eq. (49) . We 
thus have 

d 

sF[A,\u,uj\ = —5eF[A,\uj,u}\ (50) 



15 



or conversely 

5eF[A, \,uj,(j] — £sF[A,X,ou,u] 
From Eq. (50) one finds the following representation 

with the finite BRST transforms of the fundamental fields given by Eq. (49) 



{sA'^') + (sA"0 + {su'') + (so;") 



(51) 



(52) 



(53) 



sA"(a;) = gf''^^u)%x))^{x) sG)"-{x) = gf''^^uj\x)u;\x) + iaT>''\x) ■ >^ {x) 
The BRST operator s has two essential properties: it annihilates the complete action (46) 



sS'[A, A, uj^uj] =0 



(54) 



expressing the invariance of 5'[A, A, a;, a;] under the BRST transformation (49), and it's 
nilpotency 

= (55) 
Using the operator s, we now define the generating functional in the following way 



z[jj] = J: 



VA'DXDuj'Du expl -5[A, A,a;,a;] + dx 



+ • s A« + I2x • s A" + 7« + I^.sO" 



■■sX 



(56) 



Note that cu, Cu, sA and sA together with their sources J^, Jj;^, Isa, IsX are Grassmann odd, 
the remaining quantities Grassmann even. 

We proceed to vary the fields in Eq. (56) according to Eq. (49). The Jacobian of such 
a transformation is unity due to Eq. (53) as can be seen from the explicit calculation in 
Appendix A. We also know that the action is invariant under this change of variables 
S'[A, A, cu, Cu] — S[A, X, Lu', Lu']. In addition, the source terms of the BRST transformed fields 
are invariant due to the nilpotency of s and the fact that the variations are s-transforms 
themselves, e.g. S^A — esA. Only the source terms of the fundamental fields are not 
invariant and transform according to 



A«- = A'«- + d.A"- = A'«- + £sA'«- 



(57) 
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and likewise for the other fields. Thus, under the change of variables (49), the integrand in 
Eq. (56) is simply reproduced with all fields replaced by their primed counterparts and an 
additional factor 



exp < e dx 



(58) 



generated by the transformation of the fundamental source terms, Eq. (57). Because e is 
Grassmann odd we have 



expje jdx sA"-J^+sA"-J^+sa;V^+sa;V? 



1+e dx 



Inserted back into the path integral Eq. (56), the one just gives Z[J,I], which cancels the 
left-hand side of the equation. Hence, we obtain 







-I 



VADWcuVCu e dx 



exp 



(59) 



where the dots represent the exponential in Eq. (56). This has to be true for any e and 
thus the expression without e has to vanish itself. Changing the BRST transformed fields 
for functional derivatives with respect to their sources, we find the following identity 

S . S . S . S 



dx 



^^^^^^ swi 



Z[J, /] = (60) 



5P^(x) 

Finally, let us transcribe this relation in an identity for the generating functional of one- 
particle irreducible (IPI) correlation functions. To this end, we first express it by the 
generating functional of connected correlation functions PF[J, /] = lnZ[J, /]. In terms of 
W[J, I] the relation (60) reads 



Jdx 



^ , Sfe^l ™ + j,<,) ^"'1^- ^1 



(61) 



51^ ' '"^^ 5I^Jix) ' SIM ' ""'"^^ SIM J 

To define the generating functional of one-particle irreducible correlation functions, we in- 
troduce the usual expectation values for the fields in the presence of the external sources 

SW[J,I] „, , 5W[J,I] 



A^'ix) 



A"'(x) 



6Jf{x) 



SW[J,I] 
5Jf{x) 



UJ^ix) : 



5J-{x) 
5W[JJ] 



(62) 



SJ&{x) 

The minus signs in the case of the ghost fields are a consequence of our definition of the 
generating functional, Eq. (56), where we ordered the sources to the right of the fundamental 
fields. 
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Assuming that the relations (62) can be solved for the sources J, we can define the IPI 
generating functional F as the Legendre transform of W[J, I] with respect to the sources 
J. The sources of the BRST transformed fields are not Legendre transformed and play 
the role of spectators only. With the definition 



r[A,A,a;,a;;/] = JdxiA'-J^ + A"-J^" + co^J^ + a;V?] - 1^[J, /] 



(63) 



one finds 



and also 



sr 


= JX(x) 


sr 








6r 




6r 






Suj^ix) 


sr 


5W 


sr 


sw 




(5/^(x) 




SIM 


sr 


5W 


sr 


sw 




5I-l(x) 







(64) 



(65) 



which may be used to reexpress the gauge Ward identity (61) in terms of F 

sr sr sr sr sr sr sr sr 

+ 



jdx 



sr sr 

+ 



sr sr 

+ 



SA'^^ix) SI^{x) SX^^{x) SI^ Suj'^ix) S^Jx) Suj''{x) SI^^{x) 



(66) 



C. Stochastic Ward Identities 



We have included in Eq. (56) the auxiliary field A and the ghost fields uj and Co, all of 
which were not strictly necessary, but rather were included so as to facilitate our work. They 
could, in principle, be integrated out and we would be left with Eq. (11), except that we have 
now also introduced sources for the extra fields, as well as for the BRST transformed ones. 
This would suggest that there could be some sort of relations for the generating functional 
in Eq. (56) resulting from our choice to include the extra fields and sources. 

To derive these relations for Bodeker's effective theory, one starts from the generating 
functional (56), including sources of the fundamental as well as the (gauge) BRST trans- 
formed fields. Inserting the action and the BRST transforms according to Eqs. (46) - (48) 
and Eq. (53) with the definitions (6) and (41) in use, the generating functional Z[J, I] may 
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be written 

+ uj^iri^, A] + - gr'^uT^^) + A»- + cuV^ 

+ • B^'uj' + Iiyr''u;'^u;'+] | (67) 

where terms multiplying A and lj have been collected. Because the exponent is quadratic 
in the former and linear in the latter, both of these fields can be integrated. One obtains 



Z[J, I]=J DADo; 5(7') expi^Jdx 



+ ■ U^'u' + II yr^^oo'^uj"] I (68) 

with the new functionals E' and 7' defined as 

E'"[u;, A; J,, I^, 7,^] = E"[A] - iJ^ + igr'^co'r^ - aB^^'l^ (69) 
7"^[^, A; J^, 7,^] = Y[uj, A] + J? - gr'^uj'^,^ (70) 

Hence, when restricting to vanishing sources Ja—J-sA—O and J^^ — Isuj— ^ the exponent 
becomes purely quadratic in E'. Defining for brevity 

Zl[Jx, Ju!,'^s\, ^soj] = Z[Ja— 0, 3x, Juj— 0, Jo), IsA— 0, IsX, 7s^^= 0, Isijj] (71) 

we have 

Zi[J,, J^, I^, 7,^] = j m!Da; 5(7') expj-^ ^c^x E'"- E'"| (72) 

where E' and 7' both depend on A and uj as indicated in Eqs. (69) and (70). Thus, it is 
quite natural to attempt a change of variables from A and uj to E' and 7'. The Jacobian 
can be calculated in a similar manner as the Jacobian of Eq. (42), and again can be shown 
to be a constant. The resulting integral is gaussian and evaluates to a constant functional 
Zi leading to 

^i[Ja, Jw, IsX, Isuj] = const. (73) 

or likewise for Wi — In Zi 

Wi [Ja, J^, I^, 7,^] = const. (74) 
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As a consequence, any combination of functional derivatives with respect to sources chosen 
from the class {J^, Ja>, Isx, ^su/} yields zero when acting on the full generating functionals and 
evaluated for vanishing sources: 

5 5 5 



5... 5... 5.. 



W[J, I] =0 (75) 

j=i=o 



any combination 

of Ja, Ju,, Is\, IsQ 

with the same relation holding for derivatives of Z[J, I]. To obtain a corresponding identity 
for the IPX generating functional F, note that due to Eq. (74) one has on the submanifold 
defined by the vanishing of the four sources Jyi, 1^4, Juj and Isuj 

5Wa 5W^ 

A"'(a;) T T n = ^-r-h = and iJjHx) . . n = "t^A = (76) 
So r could at most depend on A, to, I^x and IgQ. However, from Eq. (64) we have 



5A''\x) ' 5a;«(a;) 

and therefore F may not depend on A or a; anymore. The same conclusion can be reached 
for the and Ig^i by looking at Eq. (65) and Eq. (75). Therefore, F must be a constant, 
this leads to 

— F[A,A,a;,a;;/] =0 (78) 

5 ... j=i=o 



^5... 5. 



any combination 
of A, u, 7sx, Isc 

which is the equivalent of the stochastic Ward identity (75) in terms of the IPI generating 
functional F. 

D. Ghost Number Conservation 

We will discuss one last symmetry of the action (46). The action is invariant under the 
global transformation 

uj"'{x) = e'"a;'"(x) 

(79) 
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of the ghost and anti-ghost fields. In addition to this, subjecting the measure "DuVoj to the 
transformation (79), i.e. to 

co%x,), u;»(x2), (I;"(x2), . . . ) (80) 
= (e^V(a;i), e-^V"(a;i), e^'^uj'^x^), e-^'^uj'^x^), . . . ) 



one finds 



(81) 



J-'(u;',a)') =det 



with the Jacobian 

/ e+^" • • • \ 

e-^" 

e+^" 

e-^" 

V ; ; ; ; ■./ 

Hence, the measure is also invariant under the transformation (79) 

TtouTiu = T>u}'T>iD' 



(82) 



(83) 



Together with the invariance of the action, this symmetry leads to ghost number conser- 
vation, which poses another restriction on the form of the generating functionals and their 
derivatives. Indeed, taking the parameter a in Eq. (79) to be infinitesimal and performing 
the corresponding change of variables 



cu^ix) = cu'^^ix) +iauj"'{x) 
iD"'{x) — [v'"-{x) — iaCd'"-{x) 
in the defining path integral (56) of the generating functional Z[J,I] yields 



(84) 



Z[J, /] = / DAVWuDu exp 



i^ia Jdx 



u;V«-u;V«+I^-sA''+PvsA"+2r„su;« 



exp 



Here we have already renamed the primed symbols again to unprimed ones after the change 
of variables has been completed. As before, the dots represent the original exponent as it 
occurs in Eq. (56). Using the fact that a is assumed to be infinitesimal, we can expand 
the first exponential and replace any fields that appear by functional derivatives acting on 



21 



the exponential (after interchanging the order of the ghost and anti-ghost field and their 
corresponding sources leading to a minus sign in either case) . The derivatives can finally be 
pulled out of the functional integral and we obtain 



Sl^^ix) 



51-^{x] 



SitSx)! 



Z[J, 7] = 
(85) 



Again, the definition W[J, I] = lnZ[J, I] implies that the same identity holds for the gener- 
ating functional PF[J, /] of connected correlation functions 



dx 



SW 



r, 



ja 



SW 



+11 



SW 



SJS ' si-^ 



7 + 1 



„■ SW 



SI: 



+ 2r 



SW 



SI"- J 







(86) 



where we have suppressed the space-time argument x and the dependence of W on the 

sources J and /. This identity in turn can easily be translated to the corresponding restric- 
tion on the IPI generating functional r[A, X,u!,Cu; I]. By means of Eqs. (62), (64) and (65) 
one finds 



dx 



ST 
Soj"- 



„ ai i ST 



Soj°- 



SI" 



SI" 



+ 2C^ =0 (87) 

This concludes our derivation of non-pcrturbativc identities for the generating functional 
(56). Explicit forms of these identities for lower N-point functions are shown in Appendix C. 



IV. DYSON-SCHWINGER EQUATIONS 

To derive the Dyson-Schwinger equations, we observe that the path integral of a 
functional derivative vanishes, i.e. 

"""ih^^^" '''' 

for any functional F[(l)]. Hence, in the case of Bodeker's theory, we obtain four different 
equations by inserting a functional derivative with respect to each of the fields A, A, cu or 
Cu into the generating functional 



Z[J, I] ^ j DA'DXDuj'Dui expj -S[A, A, a;, a)] + j dx 



(89) 



22 



A. General Dyson— Schwinger Equations 



1. Ghost (u) and Anti-ghost (u) equations 



Starting from the identity 

where the dots represent the exponent of Eq. (89), gives 



(90) 



0= / DA'DXDu'Du 



-I 

+ 9r"'{-l'U{^)-A%x) - I':,{x)-X{x) + Ilix) u'^ix) + Il{x) u%x))\ exp{(. . . )} (91) 
Expressing the fields by derivatives acting on the exponential, one obtains 



5Z 



zj»(x)=(a, + iA)-^-^r^a, 



or in terms of W 



SZ , , 5Z 



ZV-I'Uix 



5JX{x) SJ^'ix) 



SZ rh , \ 



SJ^{x) 



SJ^ix) 



(92) 



sw 



J"'(x)^(dt+^A) -^^^ - ^/'^'"^a 



5^W 



5W 5W 



sj'^{x)srj{x) 5Jt{x) sri 



X] 



-V-P^(a:) 



+ 9f 



•abc 



'rbi, ^ sw ... . 5W rb ^ ^ sw rb ^ ^ sw 



SJjix] 



Transcription to the IPI generating functional F yields 



(93) 



5T 



X] 



S'W 



8Jl{x)8ri{x, 



-Q^{x) A^\x) 



-V-i:4(^) 



•abc 



/5(x) A^\x) + /S(x) \^\x) - Il{x) u;\x) - ll^{x) u^\x) 



(94) 



The antighost equation is obtained in a similar manner and reads 

5'^W 



^^(-a. + iA)."(.) + ^/-a, 



5Jt{x)5Jj{x') 



-uj\x)A''\x') 



-gr-ll{x)u\x) 

=x 

(95) 

where x' is set to x after the space-time derivative is carried out, i.e. the derivative acts on 
the argument of Jl,{x) only. 
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2. Auxiliary field (X ) equation 



To deduce the auxiliary field equation from 







'DADXDuj'Duj 



exp{(...)} 



(96) 



wc need, among other things, the functional derivative of the action 5'[A, A, a;, a;] = 
5''°^[A, A] + S^'^°^[A,uj,uj\. However, in the present case the corresponding expression be- 
comes rather cumbersome. 

As for the Feynman rules in Appendix B, we want to use a symmetrised AA^ and AA^ 
vertex. The A dependence of the action spreads out over the three contributions to the 
dynamical action 5(°'[A, A] = ^^'"[A,A] + ^^l^^'.sIA, A] + 5;°'_4[A, A]. The corresponding 
derivatives can be written in the form 



5S'rfA,Al 



^aTX'^x) - i [5'^{adt - A) + (l - ^) didjlA'^ix) 



^^(D) _ 



(97) 



(5A'^*(x) 



= ^(-^/)V^SM^^■(^)^^^(^)^'^'(^) 



(98) 
(99) 



where V^^^^f is defined in Eq. (B.24). One then obtains in terms of the IPI generating 
functional 



5V 



5A«^(x) 



= 2aTX'\x) - i [5'^{adt - A) + (l - ^) didj\ A"^ (x) 



2! 



'^9 yabcd 
gl ijkl 



X 



6^W 



5fJ{x)5Jf{x' 
5^W 



+ A^^{x) A^'^ix') 



5fi {x) 6Jf{x) 5Jf{x) 5fi {x) 5Jf{x 



gf 



abc 



-I^{{x)u;%x)+ialUx)A'^\x) 



+ iadii:^{x) 



A'^ix) + A^\x) A'^'^ix) A'^'(x) 



(100) 
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3. Gauge field (A) equation 



Finally, coming to the gauge field equation 







DAVXDujDuj 



exp{(...)} 



(101) 



and using the derivatives 



5A"-Hx) 



[5^^{-adt - A) + (1 - ^) did,] X'^ix) (102) 

igt'" [-(1 - 1) W^dl,+ p'id, + a:)] + 2 [5^'d:+ 5'\d, + dl)] 



5A<'i{x) 

where the symmetry of Vf-^ has been exploited, together with 



(103) 
(104) 



5A"-\x) 



one arrives at 

5V 



-^r''u;\x)diUJ%x) 

K 



-i [S'^{-ad, - A) + (1 - ^) d,d,] X^^ix) - |p 



(105) 



dA''^{x) 
+2 



5Jf{x)5fi{x) 5Jf{x) 



■c/..„^ , +X''\x)A'^{x)A'^\x) 



■ A^^{x) + X'^\x) 



5Jf{x)5fi{x) ^ " 5fi{x)5Jf{x 

5^W 



[6'''d^+6'\dj + d;)] 



5f^{x)5Jf{x') 



+ X^^{x)A'%x') 



K 



S4{x') 5JZ{x) 



+ uj\x') uj\x) 



+ gf 

x'=x 



abc 



I'^{x)u\x)+i<jll{x) A-(x) 



(106) 



B. Explicit Equations for Lower N- Point Functions 

1. Definitions and General Relations 

Concerning the propagators, mixing will occur between the gauge field A and the auxiliary 
field A, resulting in four possible propagators from the gauge/auxiliary field sector that can 
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be combined into one matrix propagator. These are completed by the propagator of the 
gauge ghosts. Altogether, we define the full (connected) propagators as 

5''W[J, I] 



5J%'ix)5j'Jiy) 
6^W[J, I] 



5Jl\x)5fi{y) 



J=7=0 



J=I=0 



J=I=0 



.J=I=0 



(107) 
(108) 
(109) 
(110) 



and G'(^^)g'(x,|/) = G'(^^)5?(i/,a:) of course. In graphical representations we denote the gauge 
field by curly lines, the auxiliary field by double curly lines and the gauge ghosts by dotted 
lines. Thus, the full propagators are represented by 

X ^ y 



G^^^>^{x,y)= a,i*^QQ£iQi^^ b,j 
G'(^^)«j'(x,y)= a,i^^^Q^i2Q^ b,j 



X 



and finally 



X 



• • b 



Besides the propagators, we have to set out our definition for the self-energies. To this end, 
let us summarise the two left-hand equations of (64) in the form 

5r[A, \,uj,u;; I] 



J^x) 



;iii) 



5F«^(x) 

where the index F stands for any of the fields A or A. Taking the functional derivative of 
this equation with respect to Jaiy), where again G e {A, A}, then yields (observing that 
A, A, a; and cu are functionals of the sources J and /) 



SA'^'iz 



ST 



+ 



5^r 



X] 



SfJ{y) Suj^{z)5F'^^{x) Sfj{y) 5uj''{z) 5F''^{x) 



:il2) 



26 



Thus, using the Eqs. (62) to express the first factor in each term as a second derivative of 
W and finally setting the sources to zero leads to 



(113) 



dz 



5^r 



5^r 



J=I=0. 



j=i=o SX^'i^iz) 5F°-^{x) 

Here, the definitions (107) - (109) have been used and the terms involving ghost and anti- 
ghost fields have vanished due to ghost number conservation. 

In the following we will often encounter multiple derivatives of the generating functionals 
W and r evaluated for vanishing sources. Let us therefore introduce a shorthand notation 
where we indicate the fields with respect to which the derivatives are taken as superscripts. 
Possible LORENTZ or colour indices as well as space-time arguments appear in the order of 
the fields they belong to. For instance, we abbreviate 

5^T 



(114) 



j=/=o 



5X"'^{x) 6A^:>{y) 5u;''{z) 

In the case of VF, we also use the fields as superscripts though the derivatives are taken with 
respect to the corresponding sources, of course. 
In this new notation, Eq. (113) reads 



5^' 5'^ 5fg 5{x-y)^ jdz G^^^^jKy, z) r(^^)^,«(z, x) 



(115) 



where if is a summation index running over the fields A and A. This equation expresses the 
fact that the matrix propagator of the gauge/auxiliary field sector 

(G^^^^fAx.y) G^^'^>^{x,y)\ 



Gt^{x,y) 



is inverse to the matrix 



r?i'(x,t/) 



/r(^^)g'(x,y) r(^^)g'(x,y)^ 
\r(^^)g'(x,i/) r(^^)g'(x,y)y 



(116) 



(117) 



constructed of the second derivatives of F. Consequently, the self-energy Il'}j{x,y) is deter- 
mined via the relation 



r(^^)?(:r,y) = (A-^)(^«)«^(x, y) + n(^^)«f(x,y) 



(118) 
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where (A~^)(^'^)"j'(x, y) are the components of the inverse free propagator of perturbation 
theory (see Appendix B, Eqs. (B.6) - (B.9)), and where F,G E {X,A} as before. 
Analogously, taking the derivative with respect to of 

Sr[A,X,u,u;;I] 

and performing the same manipulations as described above leads to 

5^r 



S"'' S{x-y)^- J dz G'(^) '"'{y, z) 



bUj''{z) 8b)'' {x) 

Hence, we define the self-energy of the gauge ghosts via 



(120) 



j=/=o 



r(^")'^*(a;,7/) = - [(A-^)(")"''(a;,2/) + n(")"^(a;, 2/)] (121) 

with the free inverse propagator (A~^)'^'^)"^(x, ?/) given in Eq. (B.18). In our graphical 
representations we denote self-energies and other one-particle irreducible quantities by open 
circles. 

Though generally we are using three- vectors, in the Fourier transformation we use 
four- vector notation 

with —ikx = —ikot -|- i k • x. The proper vertex functions in momentum space are basically 
given by the Fourier transforms of the various functional derivatives of the IPX generating 
functional F. However, due to translational invariance of the theory, all these Fourier 
transforms contain a delta function expressing momentum conservation at the vertex. It 
is therefore convenient to pull these delta functions out of the definitions of the vertex 
functions. In this way, the latter become functions of one momentum variable less than 
indicated by the number of external legs. For instance, we define 

{2TTfS^{ki+k2+k^)r^''''^^'''"'j{ki,k2) = jdxdydze-'^'''-'^^y-'''^'V^^'^^^''^''^{x,y,z) (123) 

or equivalently 

Here, the two arguments of the proper vertex function r^'^'^'^^"^^(/ci, k^) refer to the (incom- 
ing) momenta along the ghost lines leaving and entering the vertex in this order. 
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The choice of the N — 1 momenta that are used as arguments of a vertex with N external 
legs is, of course, arbitrary and thereby a source of possible confusion. We therefore explicitly 
list the definitions of the other relevant vertex functions used in this work 

with k2 and denoting the incoming momenta along the G and H line respectively, and 

T^{FGHK)abcd( ,n\-- ( ^^^2 ik^{x-y)-ik3{x-z) 

i ^^^i{x,y,z,w)- J (27r)i^ (27r)^ (27r)^ 

with incoming momenta /c2, k^, k^ along the G, H and K line. Note the minus signs in 
the last two equations. The definitions above are chosen in such a way that they reduce at 
leading order to the corresponding vertices of the Feynman rules, i.e. 



p(a>a;A)a6c(^^^^^) = ^-^ P^^k^ + ... (127) 

r^^^^^tjkik2,ks) =-gV^l%k,,ks) + ... (128) 

r^'^^^^tSih,h,k,) = + ... (129) 

2. DSE for U^'^\k) 

Let us again start with the ghost equations, being much simpler than the equations for 
the gauge/auxiliary field sector. By taking the derivative of Eq. (95) with respect to u!^{y), 
one finds evaluated for vanishing sources 



(A-i)H»''(x,y) 



X = x 

J = 1 = 

(130) 

Comparing to the definition of the self-energy of the gauge ghosts in Eq. (121) then leads 
to the relation 

n('^)-(a;,,) = ^r^a,-^^^^^ 



^ ^/ 1/ I.I II 

131) 



X = x 
J = I = 



for the gauge ghost self-energy. If we carry out the functional derivative with respect to u!^{y) , 
four terms arise because any of the sources J a, Ja, and Ju, depends on cu. However, due 
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to ghost number conservation three of these terms vanish when the sources are set to zero 
and one is left with^ 



8uj\y) Soj'^iv) SJi{v) 5J^{x) 5r/{x') 



(132) 



X — X 

J = I = 



Finally, we express the connected three-point function by its IPI counterpart 

^(^c.F)a6c^^^ y^z)^ j du du'du" G^^^ "'"(u, x) G^'^^ '''''(t/, u') G^^^^f-,{z, u") T^'^^^^^''''^^ (u, u', u") 

(133) 

where F represents one of the fields A or A and G is a summation index taking these two 
values. The shorthand notation used here was introduced in Eq. (114). Note that the order 

of the ghost and anti-ghost fields in Eq. (133) is changed from W^^'^-^^ to r'^'^''^'^^ and that 
the (full) gauge ghost propagator is G'-^^"'^{x,y) = W'^'^'^^"-'^{x,y), as defined in Eq. (110). 

Now, inserting relation (133) into Eq. (132), using the property (120) of the two-point 
functions and 

yields the Dyson-Schwinger equation 

n('^) ''\x, y) = -j du'du" G^^'^^f., (x, u") ^r''' ajG^""^ '^'^'{x, u') r(--G)'=''^j; {u', y, u") (135) 

Using the definition for the momentum space proper vertex Eq.(123), we transform to mo- 
mentum space 

U^-)-\k) = - / 'If-dej^fJ G^AG)ee'fj^ _ Q^^^ ""'ik') T^^^^^^^'^i-k', k) (136) 

J [2n)^ K ■' 

The structure of the Dyson-Schwinger equation (136) is illustrated in Fig. 1. In Eq. (136) 
the field index G has a summation index taking the values G — X and G — A. In the 
graphical representation of Eq. (136) such a summation is symbolised by a solid line. This 
short-hand notation will become even more important in the other Dyson-Schwinger 
equations to follow. Thus, the right-hand side of Fig. 1 stands for two individual diagrams. 



^ It should be clear that x' is set to x only after the space-time derivative is carried out. In order to avoid 
an extensive use of brackets we decided to assume in this and similar cases some thoughtfulness on the 
part of the reader. 
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FIG. 1: DSE of the gauge ghost self-energy, Eqs. (136). Filled circles denote full propagators. 
Empty circles are used for one-particle irreducible quantities, i.e. self-energies and proper vertices. 
The solid line represents a summation of one graph with the line replaced by a gauge field and a 
second diagram with an auxiliary field instead. 



Above we have deduced the Dyson-Schwinger equation of the gauge ghost self-energy 
from the general anti-ghost equation (95). A complementary relation can be obtained from 
the ghost equation (94). By taking the derivative with respect to a)^(y) of Eq. (94), one 
obtains 



d^k' ig 



{2n] 



D 



K 



3. DSE for U(^^\k) 



We come now to the Dyson-Schwinger equations of the gauge/auxiliary field sector. 
Taking the derivative with respect to X''-^{y) of the auxiliary field equation (100) yields after 
setting the sources to zero 



5'T 



= 9rrT i)"-^ ^"-^ ^(t — ll\ yacde " 



J SX^^{y) 5Jf{x)5Jf{x'^ 



J=I=0 



(138) 



X — X 

J = I = 
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Thus, comparing to Eq. (118) one reads off tfie self-energy component 



+ [6^^dl-6^''di] 



SX'^^iy) 5Jf{x)8Jf{x') 



X — X 

J = I = Q 



acde 



5W 



3! SW>3 (y) SJf {x) 5Jf (x) (5 JJ^(x) 



(139) 



J=I=0 



To evaluate Eq. (139), we have to calculate the remaining functional derivatives and fi- 
nally transform into momentum space. Let us start with the A derivative of the connected 
two-point function. Because wc will encounter similar expressions also in the Dyson- 
SCHWINGER equations of the other self-energy components, it is useful to generalise a bit 
and do the work once and for all. Thus, with F, G and H chosen from the set {A, A}, we 
find by means of the chain rule and using ghost number conservation, together with the 
identities (64) 



5F^i{y) SJ^^{x)5J^{x') 



J=I=0 



dv 



5^r 



5^W 



J=I=0 



_5F^i{y) 5K^^{v) 5J'j^'{v) 5r^{x) SJ§{x')_ 
The field index K in this equation is summed over the two values A and A. Expressing the 
connected three-point function by its one-particle irreducible counterpart 



^,iFGH)abc^^^ 2/, -2) = ' j du du'du" ^(^^'^^^/(x, u) G^''''^%,{y, u') 

G(^^>t(.,K'')r(^'^'^>W(^^^,^.) 

and exploiting the relation (115) then leads to the identity 

5 5^W 



(140) 



j=/=o 



5F^^{y) 8J^{x)5J'^{x') 
- j du'du" G^^^'^TU^, u') G^^^'^ffix', u") r(^«'^')f;^:(y, u\ u") 



;i4i) 



Again, doubled field indices are summed over A and A (which we will assume from now on in 
all relevant cases). Finally, transforming into momentum space and inserting the definition 
of the three-point vertex function (125) yields 
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5^W 



5Fb^{y) 5J^{x)5J'^{x') 



d^k d^k' 



j=i=oJ (S^)"" (27r) 



D 



^-ik{x-y)^ik'{x-x')Q{GG'W^^j^ _ 



G^'^'^'^tUk') r(^^'^')f;^:(A;'- k, -k') (142) 
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Analogously, one can derive a general expression for the fourth functional derivative in 
Eq. (139). Using the chain rule as above, exploiting ghost number conservation and the 
identity (115), translating connected into one-particle irreducible quantities as in Eq. (140) 
and finally introducing the momentum space vertex functions (125) and (126) leads to 



(5W 



-ik{x—y) 



5E^^{y) 8Jf{x)8J^{x)5J1^{x) j^j^Q J (27r)^ ^ J (27r)^ (27r)^ 

+ G^^^'^TUk - k') G(««')f,f (^'- k") G^^"'^±,{k") r(^'«'^>;t:(^"- K -k") 

r<[L'K')h'g' y(EF' K')bc' g' /_ , /^ 
, n{GF')dc'(u _ U\ n{HG')ed' /> /_ > //x n{FH')ce' (ull^ p{L'G' H')h'd'e' full _ ,l _,in 

G'(^'^')S'(A;') V^^^'^'^%{k'- k, -k') 

, ri(HF')ec' _ U\ n(FG')cd'(y_ ,ii. r<{GH')de' yiL' G' H')h' d' e' ( ,ii _ ,/ 

+ G^^^'^fj,,{k - k'- k") G(««')f,f ()t') 



(143) 



Exploiting the identities (142) and (143) one can now readily obtain the Dyson- 
SCHWINGER equation of the II^^'^) self-energy component from Eq. (139). One finds 



1 f d'^k' 



2 J (27r) 



n i-9Wkf{^ - k', k') G(^^\Uk - k') G^^^'^ff! {k') 

jk'l'K'^ '^1 I 



1 f d k' d k" 2Yacde r,(AF')cc' (J j l\r,(AG')dd'(j l J in r,(AH')ee' (Ul^ 

2 J {2tt)^ {2tt)^ ^iklm ^ kk'V'^ tv } hi \tv tv ) Ksr mm'K'^ ) 



Y(L'G'H')h'd'e' />//_ y _yi\ (h'^ r(^F'K')bc'g' ( U _ L. _U) 



1 f d^k' d.V 



6 J {2^r (2^)^ G^^^'^'^^'ik - k'- k") G^^^'^ffik') G^^^')±,{k") 

p(\F'G'H')bc'd'e' /> / , yi_ , _u _Ul\ 



(144) 



where we have used the symmetry of the vertex V^i^ in the last three pairs of indices to 
combine the first three terms arising from Eq. (143) into one. We have illustrated Eq. (144) 



in Fig. 2. 
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FIG. 2: Dyson-Schwinger equation of the II^'^-^) self-energy component, Eq. (144). 



I DSE forU.^^^){k) 



Taking the derivative of Eq. (100) with respect to A^^{y) instead of A^-' (y) and afterwards 
setting the sources to zero leads to the Dyson-Schwinger equation for the n^'^^) self- 
energy component, namely 



(A-i)(^^)«''(a:,y) 



J=I=0 



--1 S'^' [i+adt - A) Sij + (1 - ^) didj] S{x - y) 



5 



5A^i{y) 5Jf{x)5Jf{x') 



X = x 
J = I = Q 



^9 yacde ^ 

" 3! ^klm ^^bj^y-^ sjf{x) 6Jf{x) 5JT{x) 



J=I=0 



SJf{x) 5Jf{x) 



(145) 



j=/=o 
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FIG. 3: Dyson-Schwinger equation of the H^-^^) self-energy component, Eq. (146). 



Reading off the self-energy component by comparing with Eq. (118), and using Eqs. (142)- 
(143) one arrives at 

n(A^)a.(^) ^ _ 1 (_,)^-^(k - k', k') G^^^\i,ik - k') Gi^^'^ffik') 

piAG'H')bc'd'(j^f_ 7 _7/\ 
1 r d^k' d^k" .2yacde fi(AF')cc' (J h'\ rS'^'^''>'^'^' ( h' U"-] r(^H')ee' (Uf\ 

~2 J (27r)-° (27r)-^ ^'^ >^ii' n ) ^ mm'K'^ ) 

^iL'G'H'wa;^^^y._ G(^'^')JV(^') v^^^' ^'^%C{k' - k, -k') 

1 r d^' d^k" oyacde fi(AF')cc' (J jj jjn MAG')dd' ( jn fi(AH')ee' /w/n 
6 J (27r)-° (27r)-^ ''^'^ ft ; l^ft ; mm'l'^ J 

which is depicted in Fig. 3. 

5. DSE forIl.^^^\k) 



Prom the gauge field equation (106) one obtains by taking the derivative with respect to 
A^-'(y) for vanishing sources 
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J=I=0 



s^' [{-adt - A) Sij + (1 - ^) d,d,] 6{x - y) 

-igr" [-(1 - [S'^dl^ S^'idi + 9/)] + 2 + 9/)] 

5 



Jt^A^-'Xy) 6Jf{x)SJf{x') 



2! '^^''^ 5X^i{y) 5r^{x)5Jf{x)5Jf{ 



X] 



X = x 
J = I = 



J=I=0 



-jrVlmS{x-y) 



SJf{x) 5Jf{x) 



J=I=0 



_9_ racd o " _ 

'5\^i{y) 5Mx')5J^ix) 



x'=x 
J = I = Q 



(147) 



As in Eq. (146, we have again a self-energy, a tadpole, and terms of the type in Eqs. (142)- 
(143), but we also have a new term involving gauge ghosts. It can be calculated in a similar 
way to the previous cases, and comes out to 



5Ebj{y) 5Ji{x')SJi{x) 



jDh. jDu 

" ^ " ^ ^-ik{x-y)Jk'{x-x') Muj)c'c(_y 

J {2nr (2vr)^ ^ ' 



(148) 



With this, and the previous identities, Eq. (147) can be written as 

d'^k' 



(27r)^ 
d^k' dV 



-9)V,T{^\ -k) G^'^^tk'ik - k') G'(^^')f,f (fcO 

Vi>^^^')^^,(k'-K-k') 



(27r)^ (27r)^ 



W'V^^ G^^^'^tUk - k') G^^^'^ffik'- k") G^^^-^tm'ik") 



j.{L'G'A)h'd'e' (U/_ u _un n{L'K')h' g' ( ,n ^{\F' K')hc' g' /, /_ , _u\ 



d^k' d^k" 



ig'VZif G^'^^^^Ak - k') G^^^'^t;{k'- k") G^^"'^i,{k") 



(27r)^ (27r)^ 

^iL'G'H'Wd;^^^y,_ G^^'^'^'^,^, {k') T^'^^'^'^^k' - k, -k') 

1 r d^' d^k" ig^vZiG^^^'^t'{k-k'-k'')G^^^'^f:{k')G^'^^^^^^ 

Y{XF'G'A)hc'd'e' (U, , //_ , / //n 

^ jk'Vm'K'^ ^ '^i 1 ) 



2 J (27r)^(27r)^ 



+ 



d^k' in 

% - f'^ik - k'YG^^^''\-k') G^^^^^'ik - k') ri^-^)d'c'b^k'- k, -k') 



1 r d^' . oybacd Q{AA)cdlU\ 

2 J {2ti)d ^ ''^^^ 



(149) 




FIG. 4: Dyson-Schwinger equation of the II^^'^) self-energy component, Eq. (149). 



A graphical representation of this identity can be found in Fig. 4. 



6. DSE for n^'^'^^k) 

Finally, we come to the pure gauge field component H^"^^). Because (A^^Y^^'^^j = 0, one 
has in this case 



n(^^)g'(x,t/) 



5^r 



(150) 



j=/=o 
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and thus one obtains from Eq. (106) the final identity 

n(^^)»^(A;) --J-^ {-9)V,t{^', -k) G^'^^liik - k') G(^^)ff,>') 

/jDy jDi.// 
" ^ " ^ 2T/eacd niAF')cc' ( l, _ l,>\ n{AG')dd' ( U _ n{XA)ee' ( U>\ 

(27r)-^ (27r)^ '^mifeZ feik'l'^ ft ; l^ft ft J mm'l'^ ) 

-\ J G^'^^'r^Ak - k') k") G^^-'^tik") 

^ s'Z'm'V'^ ft j U- ^j,, (^ft J i jk'r'K'^ '^1 I 

1 / g a K ■ lyeacd n{AF')cc' _ L/_ --.(AGOdrfViL/N ^(AA)ee' //^ 

2 J (27r)-^ (27r)-^ '^'^ ft J wK'^)^ mm'K'^ ) 

^{AF'G'A)hc'd'e' /. / , . //_ . _y 

+ f '1 fcda^j^ _ k'yG^-)^\-k') G^^^^^'ik - k') r(--^)'=''^'5(A;'- k, -k') 

J (27rj K 

G^'^^Zkik') (151) 

which completes our derivation of the Dyson-Schwinger equations in Bodeker's effective 
theory. 

V. DISCUSSION AND OUTLOOK 

In this work we have constructed an analytic approach to the non-perturbative physics 
encoded in Bodeker's effective theory [1, 8]. Our approach is based on Dyson-Schwinger 
equations and allows for an investigation of the non-perturbative dynamics of soft, non- 
abelian hot gauge fields that is independent of the existing lattice studies of Bodeker's 
theory [18, 19]. 

The basic starting point is to transform Bodeker's Langevin equation into a path in- 
tegral. From this path integral, in principle, one could deduce the Dyson-Schwinger 
equations. However, it would hardly be avoidable to introduce an uncontrolled gauge de- 
pendence when finally truncating these equations. To control this gauge dependence, we 
therefore enlarged the system by the introduction of gauge ghosts (which is optional in 
stochastic quantisation). This enlarged system is endowed with a BRST symmetry re- 
flecting the gauge invariance; and we have derived the corresponding Ward-Takahashi 
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FIG. 5: Dyson-Schwinger equation of the E^"^^) self-energy component, Eq. (151). 



identities. A consistent truncation of the Dyson-Schwinger equations is achieved if the 
gauge and ghost sectors are truncated in accordance with these identities. 

We also derived a second class of restrictions, so-called stochastic Ward identities known 
from stochastic quantisation [12]. These reflect the characteristic structure of the path 
integral action induced by its origin in a stochastic differential equation. 

Finally, we have deduced the Dyson-Schwinger equations of the theory. They contain, 
in principle, the possibility of (finite!) vertices coupling auxiliary fields to gauge ghosts or 
gauge field/auxiliary field vertices with more than one auxiliary field, both of which are 
not present at tree level. Whether these vertices are really non-zero, will be an interesting 
question to be decided by an implementation of our formalism. 

In combination with the gauge and stochastic Ward identities given in Eqs. (C.25) - 
(C.27), the Dyson-Schwinger equations (136), (137), (144), (146), (149) and (151) pro- 
vide all the necessary tools for an analytic study of the non-perturbative physics encoded in 
Bodeker's effective theory. In particular, it can be used to study the sphaleron rate Eq. (1), 
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where A^cs in terms of the gauge field takes the form 

Ncs{t2) - Ncsih) = dt j £x£^EtBt{x) (152) 

Restricting to the lowest correlators, wc arc then interested in the unequal time correlators 
{E^{xi)E^{x2)), {E^{xi)B^j{x2)), and {Bf{xi)B^j{x2)). The first one should approach a delta 
function, while the second one should be subleading [20]. It would be a good test of our 
ansatz if we could (roughly) reproduce the factor in front of the sphaleron rate [18, 19]. 

We close this discussion with a few comments on what is to come. Since for the hot 
sphaleron rate we are interested primarily in the infrared behaviour of the theory, the first 
thing to be done moving forward is determining the appropriate relation between the anoma- 
lous dimensions for kg and |kp. This can be done by investigating the limit when kg 0, and 
comparing with the anomalous dimension in Yang- Mills theory in three dimensions [21]. 
One can also analyse the importance of the ansatz for the vertex functions by comparing 
with time-independent stochastic quantisation [13]. 
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APPENDIX A: CALCULATION OF JACOBIANS 

Throughout this work, there appear several times Jacobians as products of change of 
variables. As is well known from the literature [15], we have claimed that they are constants 
and have generally absorbed them in the measure. To make this work more self-contained, 
we provide here a derivation of this claim. 

In order to simplify the expressions, we will suppress the colour and space indices until 
it becomes necessary. The first Jacobian that we encountered was in Eq. (5), where the 
following expression appears 




(A.l) 



40 



with 



^E[A] _ d_ 1 (5K[A] _ d_ 
5A ~di^ 2 SA ~ di 



1 + 



I f SK[A] 



(A.2) 



2 \dtj SA 

where K contains all the terms in the left-hand side of Eq.(12) without time derivatives. 
The kernel of the operator {d/dt)~^ is constrained by causality to be 0(^2 — ^i)- We then 
have 



5B\ 



Det I —j = const. • Det | 1 + - ( — 



with 



1/ d\~' 6K[A] 
2[dt SA 



ab 



(t,x;t',x') = - dt" e{t-t 



2\dtJ SA 

SKnA]{t",^) 



SA){t',x.') 



(A.3) 



(A.4) 



Since K contains no time derivatives, the functional derivative produces a delta function in 
the time variable i.e. 

sKnA]{t",^) 5^f[A](^^x) 

" - ^ ) , Abu, „/^ 



S.A^At',^') 



SA){t',^) 

where we have introduced the symbol S^ to denote a variation with respect to the x depen- 
dence only. Hence, we find 



If dy SK[(t)] 
2\di) Scj) 



(t,x;i',x'; 

Q/3 



2^^' S^A]{t',^) 



(A.6) 



Coming back to Eq. (A.3) and using Trln(. . .) = lnDet(. . .) in addition to the series ex- 
pansion of the logarithm, the determinant takes the form 



Det 



SA 



- ] — const. • exp< ^ 

^ n=l 



^ziriiixr 

n 2" 



gy^ (5K[A] 
di SA 



(A.7) 



The trace in this expression can be evaluated with the help of Eq. (A.6). One obtains 



Tr 



ay SKI A] 



dt 



SA 



dti--- dtn d"~'xi ■ ■ ■ d^~'xn e(ti - t2 
S,K^nA]{h,^2) 



4i^,"MA](t2,xO 



0(^2 - h) y '\y:yzy • • • e(t„ - 1,) "y^yyyy (a.s) 



5,X«"[A](ti,x,) 



S^AIHt^,^^) 



and thus 

, dysm 



dt 

if we set 



Jdti---dtne{tr-t2)e{t2-t3) ■■■ e{tn-ti)fn{tl,t2,...,tn) (A.9) 



f„{ti,t2,...,t„) = d^ ^xi---d 



, S^Kt^^ [A] {t2, xi) S^K^^ [A] (ts, X2) S^Kt: [A] {t,, x^ 



X. 



" S,A^^{t2,^2) S^At^its,^^) 



S.A'^i{ti,^i) 



41 



for abbreviation. Unless n — 1, however, the expression (A. 9) vanishes for any function /„. 
Therefore, only the first term of the sum in Eq. (A. 7) survives and we finally arrive at 



Det 



SE[A]\ 

— — const. • exp 
oA J 



(A.IO) 



Our next task is to calculate the functional derivative of Kf[A]. To this end, it is easiest to 
write it down in components which clarifies the structure 

+ [(l - ^) did, - 5,,a] A| + g'^r'^'^f'^^A'^^AjAl (A.ll) 

Obviously, the first term, i.e. the term quadratic in the gauge field, does not contribute to 
the functional derivative with respect to A" because it always produces a 5"'' or 5"^ that is 
contracted with the structure constants f"'^'^ in front of the square bracket. The linear term, 
on the other hand, only contributes a constant that can be absorbed into the constant in 
Eq. (A.IO). Thus, we only have to take care of the third order term which leads to 



Det 



(A.12) 

where f"-'^p'^ — CaS"-^ as usual. However, in dimensional regularisation S^~^{0) gives zero 
as a consequence of the general rules of D-dimensional integration, and the determinant is 

simply a constant. 

We came across another determinant in Eq. (34) 



Det 



(A.13) 



One finds 

5-r(i,x: 



S'^5ij5{t-t')5^-\x-^') +gf 



acd 



and thus because cu is infinitesimal 



— -[a;1C](^,x)r(t,x)] (A.14) 
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= 1+ Idtd^-^x gf 



racd 



5C'''(t',x 



—-[u^%]{t,^)C%^)\ 



t'=t 

= X 



(A.15) 



The functional derivative acting on produces a 5"'^ and therefore does not contribute 
because the Kronecker delta is contracted with the structure constants. To determine the 
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remaining functional derivative of a;'^[C], let us formally integrate Eq. (27) 

t t 
a;"(i,x) = a;«(-oo,x) - J dt" [//[A]u;]"(i", x) - J dt" Sv''[A]{t" (A.16) 

—oo —oo 

Since H[A] and t^i'^fA] are local functional in time, this equation for u; has a causal char- 
acter, i.e. u!{t, x) does only depend on the values of the gauge field A(t", x) at times t" < t. 

On the other hand, Eq. (26) leads to 

t t 
aA%t,x) = a A%-oo,^)- J dt" [ D"' x + aD"V[A] ] (t", x) + j dt" C{t'\^) (A.17) 

—CO — oo 

and A(t,x) itself only depends on the stochastic force ^(t",x) for t" < t. Hence, neither 
A(t,x) nor C(;(t,x) have a dependence on ^(i",x) unless t" < t and in taking the functional 
derivative of Eq. (A.16), we can restrict the integration range accordingly 

5a;1C](t,x) _ _ L. '^[^[A]c.]"(t--,x) _ 5 6v'^[A]{t" 



8C^{t',id) J 6C'{t',x.') J SC^{t',:>c') 



t' 



Evaluating this relation for t — t' as in Eq. (A. 15) leads to 

<5c.'"[C](f.x 



= (A.19) 

t=t' 



SC'it',^') 

The only way to escape this conclusion would be an integrand that is singular in time. 
However, if Suj/6( appearing under the integral in Eq. (A. 18) was singular, the integrated 
expression would be finite which again is 6uj/6(. Therefore, Suj/5( can not be singular. 
6A/6( on the other hand can not be singular neither because of the same argument applied 
to the functional derivative of Eq. (A.17) with respect to Thus, we conclude 

Det (^^) = 1 (A.20) 

which completes the proof. 

During the introduction of gauge ghosts to the path integral, Eq. (42), there appears in 
our work another Jacobian. We can see that it has the same form as the one we have already 
calculated, but with 

lK^[u;,A]{t,^) = -- (D»^ Vu;'){t,^) (A21) 
Hence, we can rely on our general result for the determinant, Eq. (A. 10), 

Det ^\ ' ^ = const. • exp<^ — 6 dtd^'^x — } A.22 
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The functional derivative with respect to spacial variations is given by 

^ r^ab^ -gr'^m ■ V5^(x - x') 5^"^ (A.23) 

and thus, evaluated for d — a, gives a constant because the A dependent contribution is set 
to zero due to the antisymmetry of the structure constants. Note that, this time, we did 
not have to rely on dimensional regularisation to proof the constancy of the determinant as 
we had to in the case of Det(5E[A]/5A). 

When we performed the BRST transformation in our derivation of the Ward identities, 
Eq. (49), one more type of determinant appeared. In general, if Xa are Grassmann even and 
■J?! Grassmann odd quantities, a mixed change of variables of the form 

Xa^x'^ + efaix',!^') 

^i^^', +e(f)iix','&') 
with e being a Grassmann odd parameter leads to a Jacobian 



(A.24) 



J = 1 + £ str(M) 

In this expression, the matrix M under the super trace is given by 



(A.25) 



M 




( dfa _dfa \ 



dx' 



dx' 



d-d'. 



and hence 



str(M) = tr(A) - tT{D) = ^ + ^ 



(A.26) 



(A.27) 



dx'g^ d'&[ 

(See e.g. [15], Section 1.8.2. Note, however, that in our case e is Grassmann odd which 
leads to the additional minus signs in the matrix M when e is commuted with the derivative 
d/d^). 

In our case, we have two sets of commuting variables, A"*(x) and X'^x), and two sets of 
anti-commuting ones, lj"'{x) and uj"-{x). Therefore, the Jacobian is given by 



J — 1 + e dx 



SsA" 



[X 



6 sX^^Hx) ^ S sio'^ix) 6 sLj"'(x) 



M'"*(x) ' dX'^'ix) ' Scu'^^ix) ' 510'" (x) (A.28) 
However, any of these functional derivatives vanishes as a short glance at the BRST trans- 
formed fields in Eq. (53) makes obvious: The derivative always produces a Kronecker delta 
that is to be contracted with the structure constants. Consequently, the Jacobian of the 
change of variables (49) is unity. 
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APPENDIX B: FEYNMAN RULES 



The action, as given by Eq. (46), is 



(B.l) 



with 



5(°'[A,A] = jdx (jrA"-A"-a"-(D"^x B'' + (t(A"- iD«^V-A'')) (B.2) 
5(°°)[A,a;,a;] - Jdx ^- u^u" + ^u^B'^ ■ Vcu'''^ (B.3) 



1. The Propagators 

The free, quadratic part of the dynamical action 5**°' [A, A] can be cast into the following 
symmetric form reflecting the mixing that occur between the gauge field A and the auxiliary 
field A 

with the matrix 



5r[A,A] = jdxdy ^ (A-(x),^-(x)) {A-')^;{x,y) 



(B.4) 



(B.5) 



and 



iA-')^'Xix, y) = -i [{+adt - A) 5,^ + (1 - ^) 9,9,] 5{x - y) 
-1 5-^ [{-adt - A) 5,,- + (1 - ^) - ?/) 





(A-i)(^^)^^(x,y) 



(B.6) 
(B.7) 
(B.8) 
(B.9) 



We denote by non-bold symbols combinations of time and space variables, e.g. 5{x — y) 
S{tx — ty) 5^~^{yi — y). The matrix A~^ is symmetric in the following sense 



(A-^)(^^)^^(x,y) = (A-^)(^^)5«(y, 



X] 



(B.IO) 
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Hence, the matrix propagator 

/ A(^^)^j'(a;, y) A(^^)^j'(a;, y) \ 

is given by its inverse 



\{A-^{x) X'^iy))^ (A-'ix) A'^{y))J 



jd^y A(^«)«f (x, y) (A-'r^^%(y, z) = 5^ 5., 5^" 5^(x - z) 



or equivalently 



for the momentum space functions 



FH 



g-ife(a^-2/)(A-i)(^^)g'(A;) 



(B.ll) 
(B.12) 

(B.13) 

(B.14) 
(B.15) 



(27r)« 

Note again that though we are most of the time deahng with three- vectors, in the Fourier 
transform we use four-vector notation, i.e. e~'^K^-y) — g-«feo(a;o-2/o)+ik (x-y) ig^ding to 



y-i [{+iaka+ k^) - (1 - ^) hkj\ 
In momentum space, the gauge/auxihary field propagators are given by: 



(B.16) 



a, It 



^k 



, 6, j A^^^^'}^{k) = 



% + (1 - ^) 



+i(7A;o + -|k|2 



^k 



i5' 



ah 



-iako + IkP 



+ (1 - a 
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b,j A^^^^tm 



2aTS' 



ab 



a'^kl + |k|4 



For the gauge ghosts, we have the corresponding contribution to the action, Eq. (B.3), 
comprises the free part 

S^^^'^u^u]^ y"cix(I;«(-at + iA)u;« (B.17) 

and therefore 

(A-i)(-)"^(a;,|/) = 5''\-dt + i A) 5{x - y) (B.18) 

or in momentum space 

(A-^)('^)«^(A;) = 5"^(iA;o - ^|k|2) (B.19) 
Hence the gauge ghost propagator is given by 



-k 



a • 



I) A('^)«^(A;) 



k5 



ab 



inkn — IkP 



2. The Vertices 



For the interacting part of the dynamical action (B.2) we have 

5;°' [A, A] = jdz^-igf^^X"' {l-^)A'^djA^^ + 2A^^djA'" + A''^d^A 



(B.20) 



Thus, the theory provides a 3-point vertex containing one auxihary and two gauge fields 
and a 4-point vertex of three gauge fields and one auxiliary field. To simplify explicit 

calculations, it is useful to symmetrise the vertices with respect to the two and three gauge 
fields in either case. Splitting 5'|°^[A, A] into the contributions corresponding to the 3- and 
4-point vertex 

•^SfA, A] = 5S,3[A, A] + 5S,JA, A] (B.21) 
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one obtains 



S'°l,3[^,>^hJdx^ {-ig)r''X'''[{l - ^) 5''^A'%A''^] 




where 



/■abed 

ijkl 




(B.24) 



Observing that there is an additional minus sign because we have —5*'°' [A, A] in the exponent 
of the generating functional and noting our conventions of the FOURIER transform (B.14) 
of the propagators, we find for the 3-point vertex in momentum space that is symmetrised 
with respect to the two A fields 

AA^ vertex 

a, i 



Momentum conservation is thereby to be understood. By construction, the object 
V^"^'^(k2, ka) is symmetric in the last two pairs of indices (and corresponding momenta), 
i.e. 



b,j 




c, k 



(B.25) 
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Analogously, the symmetrised 4-point vertex is found to be 
AA^ vertex 




_|_ jode jhce^^ij ^kl _ 



where F^Jf was already introduced in Eq. (B.24) and is symmetric in the last three pairs 
of indices 

j/abcd jrabdc j/acbd j/acdb jradbc jradcb fr) 

^ijkl — ^ijlk — ^ikjl — ^iklj — ^iljk — ^ilkj {D.ZO) 

For the ghost sector, the corresponding interaction term extracted from Eq. (B.3) is given 

by 

5;;;°' [A, a;, 00] = Jdx ^-zMlf-^-Cj\A' • V) a;'' = jdx tAf-bc-aj^ckg^^b ^ ^T) 

and leads to the momentum space vertex 
ujujA vertex 

a 

'■ K ^ 

h c, k 

APPENDIX C: EXPLICIT CONSEQUENCES OF IDENTITIES TO LOWER 

N-POINT FUNCTIONS 



We will now find explicit identities for the lower n-point function from the identities 
obtained in Section III. 
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1. 1-point Functions 



Let us start by explicitly writing down the consequences of Ghost number conservation, 
Eqs. (85) - (87), to the one-point functions of the theory. Taking the functional derivative 
of Eq. (86) with respect to one of the sources J^, J^,, IsA, Is\ or Isoj and evaluating for 
J = / = yields 



5W[J, I] 



SJ-{x) 
5W[J, I] 



=0 

J=I=0 
J=I=0 



SW[J, I] 



5W[J, I] 



5W[J, I] 



J=I=0 



61^{x) 



--0 



J=I=0 



(C.l) 



=0 



J=I=0 



The same relations follow for the derivatives of Z[J,I] from Eq. (85). On the other hand, 
Eq. (64) implies 



sr 



SA"-'{x) 



=0 



ST 



J=I=0 



SX^'ix) 



=0 



ST 



J=I=0 



(5a;«(x) 



=0 



ST 



J=I=0 



SO"-{x) 



--0 



(C.2) 



J=I=0 



and the combination of Eq. (65) and (C.l) gives 



ST 



Sl^ix) 



=0 

J=I=0 



ST 



Sl-Jiix) 



=0 

J=I=Q 



ST 



=0 

J=I=0 



(C.3) 



SIM 

The last first derivative of T can be computed from the stochastic Ward identities, Eq. (78) 

ST 



siu^ 

Thus, all first derivatives of T have to vanish. 



=0 



(C.4) 



j=/=o 



2. 2-point Functions 

The consequences of ghost number conservation to the second derivatives of Z[J, I\ and 
VF[J, /] are summarised in the following table, indicating for any pair of sources whether the 
corresponding second derivative (evaluated for J = / = 0) is restricted to vanish or not by 
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ghost number conservation 





Ja 


Jx 


Ju) 


Juj 


IsA 




IsLO 


Isui 


Ja 























Jx 























Jul 

























'L 























IsA 



































































































r the second derivatives of the ] 
iluated for vanishing sources J — I — 




A 


A 


u 




IsA 




Isui 


Isui 


A 























A 























CO 























CO 

























IsA 

























/sX 

























Isu> 


















































(C.5) 



0) is 



(C.6) 



In the following, we will often rely on the information summarised in these tables dropping 
certain terms that are bound to zero by ghost number conservation from our calculations 
without further notice. To start with, let us recall the gauge Ward identity in terms of the 
generating functional of connected correlation functions H^[J, /]. It was found in Eq. (61) 
to read 



dx 



(C.7) 



Taking second derivatives, a variety of possibilities arise. For instance, choosing 5/SJ^{x) 
and 5/5J^{y) yields after setting sources to zero 

S^W[J,I] 5''W[J,I] 



5Jl{y)5I2{x) 



+ 



J=I=Q 



6J%%x)6C{y) 



(C.8) 



J=I=0 
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However, due to ghost number conservation both of these terms are zero by themselves. 
Likewise, the combination of 5/5J'^{x) with 5/5J^^{y) or 5/5J^{y) does not lead to any 
new relation when ghost number conservation is taken into account. The fourth possibility 
however, combining 6/6J'^{x) and a derivative with respect to JqIv), results in the identity 



6^W[J,I] 



+ 



6^W[J,I] 







(C.9) 



J=I=0 



j=i=o 5Jt{x)5lUv) 
that will be further exploited in a moment. Considering the combinations of 5/5J^{x) with 
one of the derivatives 5/SJl\y) or 5/5J^{y) again only leads to trivial relations in view of 
ghost number conservation. The pairing of 6/6J^''{x) with 6/6J^{y) yields 



S^W[J,I] 



Ut{y)5I2{x) 



+ 



6^W[J,I] 



J=I=0 



5Jl\x)5lUy) 



= 



(C.IO) 



J=I=Q 



However, this relation is a consequence of the two simpler identities 



5''W[J, I\ 



S4{y)Sl^^{x) 







5^W[J, I] 







(C.ll) 



J=I=0 



j=i=o 5Jl\x)5lUy) 
induced by the stochastic Ward identity (75). The remaining possibilities finally, choosing 
two derivatives with respect to cu, two derivatives with respect to cu, or one with respect to 
a;, one to a) again express ghost number conservation only. Hence, up to the level of second 
derivatives Eq. (C.9) is the only restriction imposed by the gauge BRST symmetry beyond 
relations that already follow from the stochastic Ward identity or simply are a consequence 
of ghost number conservation. Imphcations of the stochastic Ward identities (75) and (78) 
are most importantly the vanishing of the auxiliary field propagator to all orders 

ir2T 



G^^^ix^y)^ '™ 







(C.12) 



J=I=0 



or, equivalently, of the (AA) self-energy component 

5^r 



5A"-^{x) 5A^^{y) 



J=I=0 



(A-i)(^^)g'(x,y) = 



(C.13) 



where in addition to Eq. (78) it was used that the {AA) component of the inverse free 
propagator is zero too (cf. Eq. (B.9)). Note that Eq. (C.13) is a special case of the gen- 
eral statement that there are no pure gauge field vertices in the theory: All proper vertex 
functions of the form 

(5"r 



r^^^-^t:ki^,y,...,z)^ 



5A"''(x) 5A^^(y) ■ ■ ■ 5A''''{z) 



(C.14) 



J=I=0 
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vanish as an immediate consequence of the stochastic Ward identity (78). Further imphca- 
tions up to second derivatives (neglecting those only expressing ghost number conservation) 
are 



5u;»{y) SI-^{x) 







5^r 



J=I=0 



together with the equivalent identities (C.ll), 



sr 



and for completeness finally 

5^r 



5W 







J=I=0 



J=I=0 



(C.15) 



(C.16) 



J=I=0 



5iUx) 5iUy) 



= 

J=7=0 



(C.17) 



j=7=o 5iU^)5iUy) 

This last identity, however, does not lead to a simple relation among the lower n-point 
functions because both of the derivatives act on sources of the BRST transformed fields. In 
general, to make sense of the above identities we will have to translate the derivatives of the 
/-type to such with respect to sources of the fundamental fields. For instance, one has 

5Z 



= j 'DAVXDujDuj {Df{x) a;^(a;)) exp|(. . . )} 



= j T)A'D\Du'Du(di(^- 

5Z 



5 



= -ft 



+ 9f 



abc 



5J-{x] 
6^Z 



9f 



abc 



s 



5J2{x) 



5 



SJ^ix 



SJ2{x)SJt{x) 



;))exp{(...)} 

(C.18) 



where the dots abbreviate the usual exponent of the generating functional as given in 
Eq. (56). Expressing this identity in terms oiW — InZ yields 



5W 
61^) 



6^W 



'5J-ix) 

Analogously, one obtains after some algebra 



+ 



6W 6W 



6J2{x)6Jt{x) 6J2{x) 5Jt{x) 



(C.19) 



5W 



512{x) 
5W 

5W 



gf 



abc 



+ 



5W SW 



:9f 



abc 



5J<^{x)5J^Sx) 5Jt{x) 5Jt{x) 

5W 8W 
5J^{x)5Jtix)^ SJ^ix) 5Jt{x) 



(C.20) 
(C.21) 



-iadi 



5W 
5Jf{x) 



+9f 



abc 



5^W 



abc 



5^W 



Uf{x) 5J^^{x] 



+ 



5W 5W 



5Ji{x) SJ^{x] 



+ 



SW SW 



SJ^{x)S4{x) SJ^{x)SJi{x 



(C.22) 
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With these substitutions Eq. (C.9) translates to 

diG^''^''\x, y) - iadjG^^^>^^{x, y) = (C.23) 
_gfacd^^{^u^A)bcd^y^ X, x) - iagf'"^W^^^^'^1]f{x, y, y) + gf^^W^^''^^'^'^{y, y, x) 

To further proceed, we express the connected three-point functions by their IPI counterparts 
and transform into momentum space. Especially note that we pull out the momentum 
conserving delta function from the definition of our proper vertices. Hence, only N — 1 
momentum variables appear in the argument of a A'^-point vertex. For instance, we use 
p(a)u)G)a6c^^^^^ ^^-^ wherc the superscript G is either the gauge field A or the auxiliary field A 
and ki and /c2 refer to the (incoming) momenta along the ghost lines leaving and entering 
the vertex in this order. Accordingly, in r^^*^^^"^J^(A;2, A;3) with F,G,H e A} the two 
arguments ^2 and refer to the incoming momenta along the G and H line respectively. 
With these definitions, Eq. (C.23) takes the form 

ik'G^'^^''\k) + ayG^^^^f^ik) = 
+ G^^^ ''\k) J gr" C^") '"''{k') G^^^^iiik - k') r^'^^^y^'fi-k', k) 

- G(^^\"^;(A;) / gMG^^^^'\k') G^'^^'^'^Xk'- k) T^^-F)d'c'a^^^j^ _ ^/ ^/^ 
J (27r)-^ L 

-ia G^^%]{k') G^^^'^f^ik'- k) r(^^^)^,;iy (A; - k', k')] (C.24) 

The indices F and G in this equation are summation indices taking the two values A and A. 
However, as we will show now, the stochastic Ward identity leads to a cancellation among 

some of the terms involved. To this end, let us express also the identities derived from the 
stochastic Ward identity in the language of full propagators and proper vertex functions. As 
mentioned above, identity (C.16) relates the normalisations of the gauge ghost and mixed 
auxiliary /gauge field propagator 

gr'^j ^ [G^^^''\k) - zaG(^^),i(^)] = (C.25) 

From the first of the Eqs. (C.15) one obtains after some relabelling 

/ T-^ af"" G^'^^^iiik') '''''{k'- k) r('^'^^)'^'«^;(A; - k', -k) = (C.26) 
J (27r)-^ 
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from the second equation 



/ 



— la 



7 9f'^G^^'^'^:'^^{k') G^^'^^ik'- k) r^^'^Yk'fik - k', k') = (C.27) 

Here we have used r(^^^'^f^^{k2, k^) = r(^"^^f^-{k3, -k2 - k^) in accordance with our defini- 
tion of the vertex functions. Let us now come back to Eq. (C.24), that was found to be the 
expression of the gauge Ward identity on the level of second derivatives. With the summa- 
tion index F taking the value A, the second integral in Eq. (C.24) consists of three terms: 
the one with the two ghost propagators and two copies of the second term corresponding 
to the two possible values G — X and G — A. The last of these terms is zero because it 
contains p*^^^^^. Moreover, the remaining two terms cancel each other due to Eq. (C.27) as 
a consequence of the stochastic Ward identity. Hence, there is only a contribution of the 
second integral in Eq. (C.24) for F — X. The first integral, however, contributes for both 
choices F = X and F = A (and hkewise if F is set to A in the second integral, G can still 
take both values G = X,A). 

The gauge BRST symmetry therefore leads to the following identity to be obeyed by the 
full propagators and proper vertex functions of the theory 

ikiG(^)ab^f^^ + ayG^^^>]{k) = 

^yacd^Hcc'^^/) G^^^^iiik - k') V^^^^^^'^'^i-k', k) 

+ G^^^ ''\k) f gr" G^"^ ""'{k') G^^'-^iiik - k') r(^-^y^'d;(-k\ k) 

J V'T^r 

- G^^^>'^,{k) I gf'''\G^^^''''{k') G^^'>''''\k'- k) r(--^)'=''-'«;(A; - k', k') 
J (27r)^ L 

-ia G^^'^^k') G(^^)J^;(A;'- k) V^^'^'^y^'r: {k - k', k') 

-iG G^^^y^^k') G^^^) j;(A;'- A;) V^^^^yi,'r\k - k', k') j (C.28) 
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